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ABSTRACT 


MoD lective O1 EMC DIOICCÎ described in this report 
was to develop computer systems which would generate the 
ement connectivity, and nodal point co-ordinates for two 
and three-dimensional finite element programs using isopara- 
metric finite elements. The computer systems and sample 


\ 
problems are discussed. 
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I. INTRODUCTION 


The increased popularity of the Finite Element analyti- 
cal method due to the availability of high speed, large 
memory computers has led to the solution of many heretofore 
Ensolvable problems in structural analysis. Although users 
of this method welcome the greater freedom given them to 
pursue more fundamental problems, they are presented the 
tedious task of generating and checking the mesh input data 
of their particular computer system. This data includes the 
arrangement of nodes within each finite element (connecti- 
vity) and the co-ordinates of all the nodes. 

Automatic mesh generation is an attempt at simplifying 
input data for the finite element programs. If the user 
can communicate with the computer in a better, faster, and 
more efficient manner, money and labor will be substantially 
saved, and the possibility of human error will be greatly 
reduced. | 

Many automatic mesh generation schemes exist, however 
these Stone ате normally propric tary "yery specialized, 
and not generally adaptable to finite elements having side 
nodes. The present system is an attempt to remedy that 
situation. 

References 3 and 6 prompted this author to devise a code 
for two and three-dimensional automatic mesh generation 


schemes, utilizing the technique of mapping from a local 





non-dimensional system of reference into the global cartesian 
Sistem of reference. 

Reference 4 is a previous work on a three-dimensional 
dutomatic mesh generation scheme using isoparametric elements, 
however the coding incorporated an interpolation scheme which 
led to ill-conditioned elements in some instances. 

Me l lustrate the effect or reduction im the amount of 
data as a result of automatic mesh generation, the following 
КООШ ре is introduced: Figure l(a) is a rather complex 
geometrical shape in the x-y plane, which can be sub-divided 
into sections (super-elements) for which one of the isopara- 
mit elements is a very close geometric approximation. 
Mese super-elements may be further discretized into smaller 
aments. The connection of the super-elements and elements 
re easily illustrated in a local non-dimensional system 
Wo ordinates, as shown in Figure 1(b), where all the super- 
Ments are squares. With the appropriate boundary descrip- 
tion of the super elements as shown in Figure 1(a), and a 
super-element connectivity diagram as shown in Figure 1(b), 
tine frid of e 1(b) can be mapped back into the 
original object O prouuce Сее ЕЕ С е Pipurer 2: 

Four automatic mesh generation schemes will be presented 
in this thesis. One is for two-dimensional linear, quadratic 
and cubic elements. The others for three-dimensional linear, 


quadratic and cubic elements respectively. 








1(a) 





Connectivity of super-elements in local system 
1(b) 


Figure 1. SEL Diagrams. 





Figure 2. Final Mesh. 
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II. 55ης И ЕЕЕ ЕСИМЕ Б 


The isoparametric concept will be discussed in this 
ember in sufficient detail to establish only the procedures 
for constructing isoparametric elements and the validity of 
Nella tionships used to obtain a mapping from the local non- 
dimensional reference into the global cartesian system of 
тететепсе. The readers are referred to references 1 and 2 
TSS further details of the isoparametric concept. 

ies serendipity family of 1soparametrieze lements I] was 
Ше теа in this paper due to the ease of transformation 
(mapping) from a local non-dimensional system of coordinates, 
Were all elements are Squares or cubes, to a cartesian 


system with curvilinear boundaries. 


TWO-DIMENSIONAL ELEMENTS 
Consider the square isoparametric elements of Figure 3. 


Figure 3(a) is a linear element consisting of four corner nodes. 





Пес υπ πα quadratic element cubic element 
3(a) 3(b) IC) 
Figure 3. Isoparametric Element Configurations 


μι ος n) Coordinates. 





Figure 3(b) is a quadratic element consis time sore IA NC 
modes and four side nodes. Figure 3(c) is a cubic element 
consisting of four corner nodes and eight side nodes: 

A curvilinear element, such as that of Figure 4, may be 
Produced by mapping the appropriate square of Figure 3 into 


Me x-y plane using the relations 
n n 
x = 2 NDX у= Е № У; (1) 


where (X. ,y;) are the cartesian coordinates of node i. 
N; are shape functions of the serendipity family, listed in 
Section 11.C.1., n is the number of nodes on the element 
er 4 for linear elements, 8 for quadratic elements, 12 
for cubic elements). 

Note that the convention for node numbering within any 
ment (See Figures 3 and 4) 1s to begin at (€,n) = (1,1) 


and to number the nodes proceeding counterclockwise. 





Figure 4. Curvilinear Element Mapped 
in Cartesian Coordinates. 
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Er THREE-DIMENSIONAL ELEMENTS 

Linear elements, Figure 5, have eight corner nodes. 
Quadratic elements, Figure 6, have eight corner nodles Arnd 
meme side nodes. Cubic elements, Figüure nave eight 
Seamer nodes and twenty-four side nodes. 

Similarly, as with two-dimensional elements, three- 
dimensional curvilinear elements may be produced (mapped) 


using the relations 
\ 


where (X: ,Y:521) are the cartesian coordinates of node 1. 
Ecuwbown in Figures 5, б апа 7, the convention for node 
Exubering within any element is to begin at (£€,n,2)= (1,1,1) 
and to number the nodes proceeding counterclockwise about 


Per zeta axis. 
6(-1,1,-1) 


БТЕ 


Ш (0, -1,-1) 





4 ΠΠ] 


Figure 5. Three-dimensional Linear Element. 


11 








Figure 6. Three-dimensional Quadratic Element. 





Dunne mmecsdimensdonal Cubic Element: 


e 


IN? 





C SHAPE FUNCTIONS 


Bet. Eo = Eur yc s qup S 


O 1 O 1 66; 


1. Two-Dimensional Elements 
a. Linear elements 
4 


осе Еш 5(4) гот identirication orf the nodes: 


N; = (1/4) (E, +1) (n, +1) (3) 


where 54 = +] and ΠΠ; El, 


b. Quadratic elements 


See Figure 3(b) for identification of the nodes. 


Ecrger nodes: 1,3,5,7 
N. 7 (1/4) (1*£,) (1*n,) (£,*n,-1) (4) 
where 5i = +1 апі п; = tlhe 


pede modes: 2,4,6,8 


=, 
1 


(1/2) (1-5) (1 +п,), (5) 
for nodes 2 and 6, where > +] 


2 


М = (1/2) (1+8) (1-1), (6) 


for nodes 4 and 8, where a e 


c. Cubic elements 
SSI (CO de Ear AO nodes: 


Toner nodes: 1,4,7,10 


N; = (1/32) (1+ €,) (1 +n.) (9 (E%#n*) - 10] (7) 


where 5i - t], ΠΗ +i] 


15 





Edo nodes: 5,6,11,12 
Z 7 
N, = (ϑ/52}{1:6.}{1-η΄)(119η ), 


where 5i τε]; 1 5 аА. 


pide nodes: 2,5,8, 9 


N, = (9/32) (1+n ) (1-65) (1+98,), 


where $8 = +1/3, υπ... S: 
\ 


ph. Three-dimensional Elements 


л. Linear clements 


See Figure 5 for identification of nodes. 


Ny = (1/8) (145) (l+n у, 


b. Quadratic elements 


Cee Ficurée © for identilic cation of nodes. 


Μι Ποάξε5π: 1,3,5,7,13,15,17,19 


1 


where - = +1, D m Em qu 


Side nodes: 2,6,14,18 


N, = (1/4) (1-65) (1*n,) (1+2,) 


where 5i Ξ 0, ams F IU а SENT 


Side nodes: 4,8,16,20 


М, = (1/4) (1-12) (1+) (1+5) 


14 


(8) 


(9) 


(10) 


ЕО ο Саа ΠΡ 


(12) 


(13) 





where Si = Xll nc πα n ss 
fide nodes: 9,10, 11.12 


N, 7 0/4) si. (1*n,) (14) 


where - = +1, йт + ο. 0 


c. Cubic elements 


See Figure 7 for identification of nodes. 


Rer nodes: 1,4,7,10,21,24,27, 30 


N, = (1/64) (14E,) (14ng) (+65) D9 (E^ n^ z^) -19] 
(15) 


where S; ΕΕ r] ες. = +] 


e nodes: 2,3,8,9,22,23,28,29 


М; = (9/64) (1-5) (1+9Е ) (1+1) (1+6) (16) 


where ος М5, Da, ce 


ο 5,6,11,12,25,26,31,32 


N; = (9/64) (1-n7) (1+9n,) (1+) (1+5) (17) 


where E; Fi |D = Аа 


1 
Side nodes: 13,14,15,16,17,18,19,20 
N; = (9/64) (1-2) (1#9¢,) (14E,) (lng) (18) 


1 


where E = Jg E tally aor, 
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III. BASIS OF AUTOMATT Gale she GENERATION SCHEME 


E mentioned зп {πε introduietionmmrc terencer 5 prompted 
EE utilization of the following scheme in this thesis 
research. 

The essence of the scheme is the use of simple mapping 
nections to transform simple orthogonal connectivity dia- 
mM: into real geometrical shapes. Considering the partic- 
ular two-dimensional case of a parabolic quadrilateral of 
Eu Ure A, in which the x and y co-ordinates of eight nodes 
Encoknown, we can write equation (1) from Chapter II, where 
N; K Shape function associated with each node and defined 
ИП terms of a local non-dimensional coordinate system € and n 
which has values ranging from 1 to -1 on opposite sides. 

Typical shape functions are given in II.C., and for this 
ple equations (4), (5), and (6) apply. If the co-ordinates 
E TE nodal points are known, then cartesian co-ordinates of 
any specified point (&,n) can be found by equation (1). 

If the whole ЕРА in which ΗΠ to be generated 
could be described adequately by a quadrilateral of the shape 
in Figure 4, a mesh of any refinement could be automatically 
memetated inside it by specifying the coordinates of the 
eight boundary nodal points and the number of sub-divisions 
IM the E and n directions. The above scheme is shown pic- 


torially in Figures 8(a) and 8(b). 
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BEmodratic Super-Element 
local coordinate 


Quadratic ирст Еее 
Arte ones ood mae 


system system 
8 (а) 8 (b) 
=o 
D S Point ot inflection 
εν ο. E "Mir Parabola 
ТС 
== 9 





Super-Element Selection 


ο αν 
EMper-Element уре | 
A quadratic 
B linear 
C linear 
D cubre 
E quadratic 


corresponds to corner nodes 


corresponds to side nodes 


col 


Боли ерои е саше nodes for a single Super-Element 
whose adjacent Super-Element is of a different type 


Figures 
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If the whole region could not be adequately described by 


mroure 4, then we could subdivide the region into more than 


Ene quadrilateral and apply the same principles to each as 


meecified above. These region 
Bo as Super-Elements hereafter 

Mie local Super-Element ot 
Meastically distorted, even up 
lie along the same line in the 
triangle or semi-circle. Care 
make any corner angles greater 


may result [1]. 


ced between corner nodes tO 


subdivisions will be referred 
(see Figure 8). 

Ficuremo(a) can be quite 

to a point of making two sides 
cartesian system such as in a 
must however be taken not to 


than 180° as a non-uniqueness 


In addition, side nodes should be equally 


insure equal spacing of sub- 


divisions in the & and n directions and thus the x and y 


directions in the mapped cartesian system. 


The criteria above also apply to the linear and cubic 


isoparametric Super-Elements for two-dimensional meshes, and 


E near, quadratic and cubic isoparametric Super-Elements 


for three-dimensional meshes. 


The decision on whether to 


uSéE теат qguadratic or 


cubic Super-Elements will depend entirely upon the region 


boundary to be mapped. 


Linear Super-Elements will map 


quadrilaterals with straight sides; quadratic Super-Elements 


КИШ map quadrilaterals with straight or parabolic sides; 


ANS ouUper-Elements will map quadrilaterals with straight, 


Parabolic or cubic sides. 


Figure 8(C) 1S an Illustration ot 


EMper-Element selectrom based upon region (structure) boun- 


dary conditions. 
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The final generated mesh for 3 given Structure may Con: 
tain nodes on the boundary which do not exactly coincide 
with the actual structure boundary, or the elements within 
May not be represented sufficiently accurately. In this 
case it may be necessary to have some of the coordinates 
majusted once the mesh has been generated. This adjustment 
will be quite easily accomplished with the aid of computer 
mots, which can be obtained as output to the mesh generat- 


ing programs discussed in Chapter IV. 


ES 





IV. DISCUSSION OT COT UIT TRONS 


ihe computer programs presented in this thesis were 
του το support the computer systems 'PLISOP' and 'TRISOP', 
both of which were coded by Professor G. Cantin, and his 
meudents of the Naval Postgraduate School, Monterey, Cali- 
πια. TRISOP performs a structural analysis of three- 
dimensional problems using isoparametric finite elements, 
mite 'PLISOP' performs the same analysis of two-dimensional 
problems. The input data to each system is quite extensive 
with regard to element connectivity and nodal point coordi- 
nates. 
| The automatic mesh generating systems discussed in this 
chapter compute the element connectivity and nodal point 
we dinates required for the above two computer systems. 
Four mesh generating programs will be discussed. The first 
Ki rates data for 'PLISOP' with linear, quadratic and cubic 
elements, and is called 'PLIMEG'. The other three programs 
Benerate data for ' TRISOP' with incar mana οἵπο απο ИВЕ 
elements respectively. They are called 'TRIMEG-1', 'TRIMEG-2' 


Ende! TRIMEG-3'. 


ΠΕ GENERAL PROGRAM 
Wile α simple mesh generator using only one super. 
element may be sufficient for some purposes, it is limited to 


Simple geometries with but a single uniform property. For 


` 
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memplete cenerality, as suggested at the end of Chapter TII 


the scheme was extended to many super-elements, either linear, 


τας, cubic or any combination of the three, and repre- 


K d in a 'checkerboard' pattern, where each of the Super 


Elements may be sub-divided and each may define a material 


a single property. The above scheme is presented in 


sections B and C of this chapter for two and three-dimensional 


Mees respectively. 


B. TWO-DIMENSIONAL MESH GENERATION (PLIMEG) 


l. 


8. 


Definitions (( see Figure 9) 


ОЕШ eal acronym, wsed со refer to a super: 
element as defined previously on page 7. 


Mas fer Grid: A 10 x 10 G€heeker-board array 
of SELs which is utilized as a working 
board for SEL numbering and orientation. 
The master grid SELs are numbered sequen- 
Оа у along grid COLUMNS 


Rows: me ubdivisions of a SEL in nñ direction 


ο πας i subdivisions of aSEL INM = direction 


Elements: those- portions ofr oa ΕΙΡ τευ Dy 
rows and columns 


per element Selection 


The number of SELs required in a given mesh is deter- 


mined from boundary considerations and the types of materials 


from which the object is composed. Enough SELs must be 


incorporated to adequately define the boundary of the object, 


wat least one SEL must be utilized for each different type 


ut material. 


Object. 


Figure 10 shows an example for a particular 


Tes ασ πμ note in Figure 10a) thatthe 


common boundary of SEL 1 and SEL 2 is based on material type, 
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SEL #12 ο που DOIN) 
Master Grid 9 (b) 


9 (a) 
A indicates corner nodes; (jindicates side nodes 


All nodes indicated comprise SEL boundary nodes- 


Figure 9. PLIMEG Definitions. 


while the remaining boundaries are determined from geometri- 
cal considerations only. 

mec [he Ποπ апа type Of SELLS Mave been deter- 
mined, the proper SEL orientation in the local nondimensional 
system of reference must be decided upon, as shown in Figure 
NED). The local system 1S now superimposed upon the master 
grid of Figure 9 to obtain the SEL numbering scheme of 
Figure 10(b). Note: this superimposition may be anywhere 
within the boundary of the master-grid. 

oe Mesh Criteria 

'' ΠΠ ΠΗ: musi πε ΠΡ ΕΙΠΕ". 

a. All elements within a given mesh must be of the 
Sue type (i.e., linear, or quadratic or cubic), and only 


mice type Of SEL may vary. 
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Material 


Material 
j. 


Cartesian SEL Diagram 


10(a) 


Corresponds to corner nodes 


corresponds to side nodes 





corresponds to side nodes for 
sample SEL whose adjacent 


BEL is of a different type Local SEL Diagram 


TO(E) 
1 To Columns (linear) 
EN rows, 5 columns (quadratic) 
WES τους, 7 columns (quadratic) 
uc 5 rows, 6 columns (cubic) 





Те И 
TOE) 


UG TCU ee LIMEG Example. 
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b. The number of rows in all) SELS of any master grid 
ES must be equal τα, SELS 2 ο ο. 5. 

с. The number of colima Eso ia Cen 
Beer grid column must be equal ι ο SELs 11,12,13...20). 

d. SELs may be connected to adjacent SELs in the 
following ways (see Figure 11): totally connected, connected 
at one corner only, totally disconnected, and connected at 
Mier corner nodes only but not along common edge. Figure il 
shows only horizontally connected SELs, but the same connec- 
meron schemes apply to vertically connected SELs as well. 

4. Bandwidth, Node and Element Numbering 

Elememtznumberıng 15 always alone master grid 

Columns, and with the exception M discontinuities between 


cent SELs, such as in Figures 11(b) and 11(d), node 


O | 
Ab 





4 


Totally Partially Totally, 
Connected Connected Disconnected 
EE 11 (0) E) 


ο Έτος aC 
СОСТ Ое И 
τ." τα} 


ντε, ΠΕΡ Connections. 





mUDbering is-also along master prid columnis E D v Io, EI. 
obtain a node numbering scheme which results in the best 
Ap tational etfticiéncy, sineew PETSOP and Бор Far 
designed to be used with a banded stiffness matrix solution, 
mers prudent that the number of rows in the master grid row 
memcaining the most columns be less than or equal to the 
total number of columns in that particular master grid row. 
Mhe reader is rtg to Section O bar hiii P chapteY вот 
additional information regarding this topic.) For example, 
Ш роте 10 the master grid row containing the most columns 
EN ud row two, in which there are sixteen columns, while 
ἵπς number of rows in this grid row is five, thereby satis- 
fying the above criteria for bandwidth. Another example 
Boncerning bandwidth is as follows: suppose the object to 
E Considered 1S in the shape of an "L" and that all SELs 
Seltain an equal number of rows and columns (1.e., 3). Two 
hl rent SEL combinations could be utilized to describe 
Me mesh as shown in Figures 12(a) and 12(b). However it is 
Æ dent from Figures 12(c) and 12(d) that bandwidth would be 
Fi ater for the shape of Figure ο. ος στο are six 
rows in master grid column one for Figure 12(d) vice three for 
ure 12(c). 

There will be many instances when the geometries of 
Figure 12 will present itself, and when this condition arises, 
tne user must determine whether element shape within the mesh 


cag mesh with minimum bandwidth is the more crucial. 
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Cartesian SEL Diagram Cartesian SEL Diagram 


12 (а) 12 (b) 


Pu — - res nn mn w 


kocal SEL Diagram 
ПА Се) 





Local SEL Diagram 
12:085) 


τς το. (The "p'un papa 


Whenever discontinuities such as those of Figures 
11(b) and 11(d) exist in a mesh, a good guideline to follow 
lS to have them exist above or below a SEL, if possible, 
Eur than to the left or right, since the latter type of 
lV Connection occasionally results in a somewhat non- 
efficient node numbering scheme due to the complicated 
ιτ incorporated in Subroutine CONN of PLIMEG. 

EC PEIMEG Output ) 

πο major portion of thefonmtput witi PLIMEG ας 
ment connectivity and nodal coordinates, however addi- 
tional output has been incorporated to aid in the use of 
PLISOP. To begin with, all data input to PLIMEG will be 
meproduced as Output as a check to the user, in event a 


mispunched data card is experienced. The total number of 
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di ments and joints in the mesh together withthe half 
Eundwidth will also be printed. -In addition, an option to 
Munch on data cards the element connectivity and coordinates 
Sf the joints is incorporated. Lastly, a plot of the mesh 
may be obtained on the printer, whereby Subroutine GRID calls 
Enbroutine UTPLOT in the IBM 360 source library. However, 
расе plots on the printer are somewhat distorted and only 
wn ast ot points, a plot can also be obtained on the off- 
line plotter (Calcomp) whereby Subroutine GRID calls Sub- 
momtane DRAW in the source library. The resulting plot is 
(jose shown in Figures 2 and 10(c). The printer plot 
ШОШ be used in determining proper mesh discretization and 
then Calcomp for the final plot of the mesh. Both plots can 
ШОШ Бе obtained on the same computer run. 

Appendix A discusses input data preparation for 
PLIMEG while Appendix B discusses sample problems and shows 
the input/output formats, the node and element numbering 
πο, and the resulting plots. These two Appendices may 
be removed and utilized as an instruction manual for PLIMEG 
once the contents of this and the previous chapters are 
comprehended by ESAS e me 
INS TAREE-DIMENSTONAL MESH GENERATION (TRIMEG) 

IRTMEG "actually consists ot three separate programs: one 
mr each type of element (i.e., linear, quadratic, or cubic). 
sS are called TRIMEG-1, TRIMEG-2, and TRIMEG-3 respectively, 
However, since the computer scheme and the data input incor- 


perated in these programs is identical, they will be discussed 
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stone (1.g., TRIMEG),  TRIMEG also incor orates many Of 

the ideas of PLIMEG, therefore only the dissimilarities 

between the two schemes will be presented in this section. 
IW Definitions (see Figure I3) 

a. Master Grid: a 5x5x5 three-dimensional "checker- 
mena array of SELS. The master grid SELs are numbered 
sequentially along master grid columns in successive master 
mma Slices. 

ν᾽. Trees mer dıvısıons of SEL in geda rection 


ОРЕ еШ Пы nose portions of a SELereated Dy 
rows, columns and slices 


Mesh Criteria 

Unono Mowing rules must be adhered to at all times. 

a. The number of rows in all SELs of any master grid 
ior MUSE be equal ڪ,.€.1(‎ 1,6...21,26,31...46,51,56...71,76,81, 
КО ОТОТ, 106, ...121). 

SS Шс лиш ет ο К кш iņ all SELS of any master 
s column must be equal (1.e., 1-5,26-30,51-55,76-80, 
η] - 105). 

c. The number of slices in all SELs of any master 
EE -:ce must be equal (1.e., 26,27,28...50). 

ο πο ασ ο adjacent SELSmayZberconnected In 
any of the ways shown in Figure 11 (applies to all six faces), 
ce exception of Figure 11(d). If this situation exists 
in a mesh, then four SELs must be utilized as shown in 


Iure 14. 
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Figure 14. TRIMEG SEL Connection fOr 
Similar Connection o£ 
Figure ΠΠ 


3. Bandwidth, Node and Element Numbering 

Node and element numbering is identical to that of 
PLIMEG for each successive slice in the mesh. The following 
rules apply to obtain minimum bandwidth. The number of rows 
in the master grid row containing the most columns should be 
less than or equal to the total number of columns in that 
particular master grid row. In addition, the number of 
columns in this master grid row should be less than or equal 
ШОШ е total number of slices in this particular grid row. 
ШЕЕ 5 спете сап be visualized by sketching a single SEL with 
one row, two columns a three slices and then numbering the 
nodes along columns in successive slices. The half-band 
width is 30 as shown in Figure 15(a). Using two rows, one 
column and three slices, produces α πο ποπ TOE 3355 
ENS The reader should sketch several possibilities for 
mew, column and slice combination (Figure 15) to convince 
Iummself that orientation of the mesh determines bandwidth, 
and thus efficiency, when solving problems with PLISOP and 


TRISOP. As can be seen in Figure 15, numbering is identical 
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|x 2x 1; 
half-bandwidth = 30 half-bandwidth = 33 
15(a) 15(b) 





1η 
ον αὶ XIS 
half-bandwidth = 54 half-bandwidth = 36 
USES Ec) 


ΤΕΙ. 
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Eo that of PLIMEG for a single SEL withioncssl i CDS ОЕ 
the third dimension. For PLIMEG the half bandwidth for 
cure 15(a) would be 8. 
4. TRIMEG Output 

TRIMEG output is identical to ELIMEG except for the 
plotting package. With TRIMEG a plot can only be obtained 
Ec calcomp plotter. The plots are of three different 
MOojJections on the planes produced by rotating the mesh 
about the Z-axis, the latest x-axis and the latest y-axis 
aeo C tively. The three angles of rotation (Euler angles) 
and the resulting axes are shown in Figure 16 and the scheme 


is contained in subroutine TRFR of TRIMEG. 





апре οὐ »rotationmabout 22 axis 


= angle of rotation about x, -axis 


τ, Ὁ 
| 


II 


angle of rotation about y, axis 


Um WE xesSRotatrions for Projection Plotting. 


о 





Appendix C discusses Input data рер за Кето or 
TRIMEG while Appendix D discusses sample problems and shows 
the input/output formats, the node and element numbering 
Schemes, and the resulting plots. These two Appendices may 
be removed and utilized as an instruction manual for TRIMEG 
Eu Che contents of this and the previous chapters are com- 


Drehended by the user. 
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V. CONCLUSIONS 


e utilization of PLIMEG and TRIMEG automatic mesh 
generation schemes along with existing finite element pro- 
grams, now makes 1t much more practical to solve realistic, 
complicated problems, since the amount of input data has 
meen reduced as much as 95 per cent. This decrease in the 
ШШШ ОГ input data cards required for the finite element 
programs not only reduces the chance for error, but makes 


momen SOlving more of a challenge and not just a headache. 
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VI. RECOMMENDATIONS 


IO OF STRUCTURE MESH 

inaddition to the present plotting routine incorporated 
Il IMEG, it would be desirable to obtain a plot of the 
5165} which did not contain hidden lines. Profes- 
sor ште» Cantin of the Naval Postgraduate School, Monterey, 
California is presently working on such a solution. Adding 


Meee tual node numbers to the computer produced plots would 


also be of great value. 
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APPENDIX 8 


PLIMEG EXAMPLES, INPUT AND OUTPUT 


This appendix discusses two example problems. The first 
as how SEL selection and connection procedures dealing 
with an eight SEL object with discontinuities along some 
edges, and also the node and element numbering schemes incor- 
porated in PLIMEG. The second will be a more straightforward 
example illustrating the detail of input data and output 


formats. 


A. EXAMPLE ONE (Padeye; see Figures 17 and 18) 

The common boundary of SEL 21 and SEL 22 is based on 
material type, while the remaining SEL boundaries are deter- 
mined from geometrical considerations only. As noted in 
eure 17/(a), the common boundary of SEL 2 and SEL 42 must 
be connected manually once the mesh has been generated. 

This is accomplished as follows (see Figure 17(c)). The 

БТ еер of Figure 18 will not contain nodes 77-80, there- 
fore remove these nodes from the output punched deck of 
coordinates and change the element connectivity cards for 
Eenents 55,56 and 57. In other words, node 77 becomes 

node 1, node 78 becomes node 2, node 79 becomes node 3 and 
node 80 becomes node 4. The number of nodal points (joints) 
will now be 76 vice 80, while the number of elements will 
Γη the same (1,e., 56). Half bandwidth will be increased 
from 20 to 148, but may be held to 50 if SELs 1,2,3,4,5,12,13 


and 14, for example, are used to describe the mesh. 
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Mater als, 






Material 1 


Connected 


17(a) 





2 16 18 
Local SEL Diagram Showing SEL Numbering and Corner Nodes 
17(Ф) 
0 Li ο ο д D 6] 





8 16 31 38 45 = js {51 


Subdivided local SEL Diagram Showing Node and Element 
Numbering (linear element) 


NES) 


ΤΙ 
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Бы рат сш ү lita Mesh for Example One. 


The remaining variables for input to PLIMEG will be shown 


in the more straightforward example which follows. 


B. EXAMPLE TWO (See Figures 19, 20 and 21) 

Figure 19 shows a simplified mesh generation scheme con- 
Sisting of two SELs, six elements (linear) and 12 nodal 
points (joints). Although SELs connected vertically could 
have been utilized, this would have established a mesh with 
three rows and two columns resulting in a half-bandwidth 
greater than that of the existing mesh. 

Although the mesh is quite coarse, it adequately describes 
Ehe data input for the given boundary situation. The SELs 
mould have been subdivided to produce a much finer mesh, 
however the amount of input data would remain unchanged. 

Figure 20 shows in detail the input and input format 
mecuired for PLIMEG as described in Appendix A, while Figures 


πα 22 present the output and output format. The last 
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s columns of the connectivity matrix (Figure 22) arc BH I ID 
КИШ РЕ. These quantities represent the kind of clenent 
m (he type of material for that element, respectively. 
KIND equal to one implies linear elements, two implies 

iia ratic, and three implies cubic. This information is 


npewmeired input for PLISOP. 
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(2,2) & i _ = | 2 =>, o 2) 
Cartesian SEL Diagram 
19(a) 

1 о Αἱ 5 


SET A Oy 
P10 SEL 11 - Quadratic 






E TE 3 
Local SEL Diagram 
19(b) 





18(b) Subdivided ТЕ 2 =. же: 
(linear elements) MORTE Em 
19(c) 19(d) 
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E 0 78901235456789012345678901254567890123456789012345607890 


MINE G* EXAMPLE TWO** 


J.R. ADAMEK 
2 4 0 0 
1 2 Z 1 2 4 3 1 4 
11 2 1 5 4 6 5 1 8 
1 3 20 2 
1 1 2-0 20 
1 2 4.0 20 
1 4 3.4 5.2 
11 5 2. 6 ο. 
11 8 4 
11 3 Е 807 
11 i 4 
11 4 3.4 δ. 
11 9 4 
11 6 3.4 5.6 
11 10 206 ο 0 
3.0 2 


i 1 50769012345678901234567890123456789012345678901234567890 


Iuyvure 20. PLIMEG Input Data (example two). 
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APPENDIX D 


МЕЕ T EXAMPLES SI PUT GAN DSO UM Pita 


Since the concept for TRIMEG is very similar to that of 
PLIMEG, only the figures illustrating two examples will be 
sewn in this appendix. The user should be completely 
familiar with the contents of Appendix B (PLIMEG Examples, 
ΠΤΙ ana Output) prior to utilizing this appendix and 


Appendix C (TRIMEG Input Data Preparation). 


DENEN AMPLE ONE (Crack Propagation Problem) 

Figure 23 shows the local SEL diagram indicating SEL 
numbers and SEL corner node mE x The cartesian SEL 
diagram is not shown since it would be identical to that of 
the local diagram since the object is a cube. Figure 24 
shows the final mesh, where each SEL contains 2 rows, 2 
columns and 2 slices. 

Example two will completely illustrate all input, output 


“respective diagrams. 


- B. EXAMPLE TWO (Folded Plate) 

Figures 25, 26 and 27 show a simplified mesh generation 
scheme consisting of three SELs, 27 elements (linear) and 
eighty nodal points (joints). Although the mesh is quite 
coarse, it adequately describes the data input for the given 


boundary situation. 
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Figure 28 shows in detail the input and input format 
meauired for TRIMEG as described in Appendix C, while 
ieee s 29 and 30 present the output and output format: 

Each SEL contains one row, three columns and three 


slices. 
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Flnal Mesh 


νυ Folded Plate. 
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TRIMEG Output (Input Data Check). 


Figure 29. 
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